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Abstract 

In this paper we present a convergence analysis for the Nystrom method proposed in [Jour. 
Comput. Phys. 169 pp. 2921-2934, 2001] for the solution of the combined boundary integral 
equation formulations of sound-soft acoustic scattering problems in three-dimensional space. 
This fast and efficient scheme combines FFT techniques and a polar change of variables that 
cancels out the kernel singularity. We establish the stability of the algorithms in the norm 
and we derive convergence estimates in both the and L°° norms. In particular, our analysis 
establishes theoretically the previously observed super-algebraic convergence of the method in 
cases in which the right-hand side is smooth. 

Key words: Acoustic Scattering, Boundary Integral Equations, Nystrom methods, FFT MSG: 
65N38, 65N35, 65T40 

1 Introduction 

In this paper we present a full convergence analysis of the high-order Nystrom method introduced 
by Bruno and Kunyansky in [5l [6] for boundary integral equations (BIE) related to the scattering 
of time-harmonic acoustic waves in three dimensions. In particular, our analysis establishes theo- 
retically the previously observed super-algebraic convergence of the method for smooth right-hand 
sides. To the best of the authors' knowledge, this proof constitutes one of the few instances of 
analysis of a Nystrom type method for a boundary integral equation in three dimensions. 

As is well known, Galerkin methods for BIE of the first kind have enjoyed thorough theoretical 
analyses since their inception — on the basis of ellipticity properties and discrete Fredholm the- 
ory. Compactness arguments can also be used to establish convergence of Galerkin methods for 
equations of the second kind. Few results exist, on the other hand, concerning convergence for 
three-dimensional BIE collocation methods — in which finite-element bases are used for approxima- 
tion, but testing relies on point sampling. This state of affairs has led to the widespread perception 
that, being even more "discrete" than collocation schemes, Nystrom methods for BIEs of the sec- 
ond kind with weakly singular kernels could not be easily analyzed. This paper will hopefully help 
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dispell this view and lend additional credibility to Nystrom methods — whose qualities can be very 
attractive in practice [5j. 

We thus consider the numerical method [5] for the second-kind combined field integral equation 
associated to the problem of sound-soft scattering by a three dimensional obstacle with smooth 
boundary S. The method relies on a series of geometric constructions: 1) Representation of the 
surface by means of a set of J overlapping smooth charts (parametrizations); 2) A smooth partition 
of unity subordinated to these charts which decomposes both the overall integral operators as well 
as the solution of the BIE as the sum of contributions defined on the parametrized patches; 3) A 
family of floating smooth cut-off functions rjs : S x S (see ()2.10p below) that is used to isolate 
the singularity of the kernel function, and thus produces a splitting of the integral operator of the 
form K = -fCreg,(5 + -fusing, (5 — whereby the regular operator is an integral operator with a smooth 
kernel, and the singular part which enjoys a reduced domain of integration but contains the weak 
singularity of the integral kernel. 

The Nystrom method under consideration is based on approximation of the integral operator 
by a quadrature rule which treats the regular and singular parts i^reg,<5 and i^sing,<5 separately. The 
J local charts mentioned above are used to push forward J uniform grids from the unit square 
to the surface; these grids are used for both, approximation and integration. The quadrature rule 
used for the regular operator Kj-^g^s is based on application of the two-dimensional trapezoidal rule 
in the parameter space [0, 1] x [0, 1] for each one of the J parametric domains: since the corre- 
sponding integrands are smooth with compact support strictly contained in the unit square, these 
trapezoidal quadratures give rise to super-algebraically accurate approximations. For the singular 
part, a change of variables to polar coordinates around each integration point is applied, which re- 
sults in a smooth integrand to which, upon necessary Fourier-based interpolations, the trapezoidal 
rule is applied for radial and angular integration — once again yielding super-algebraically accurate 
approximations. 

As a result of these constructions we obtain an algorithm that evaluates the action of an 
approximating discrete operator on a continuous function, using only the values of the function 
on the selected quadrature points. The algorithm can then be accelerated by some of the available 
acceleration techniques: equivalent sources and FFT techniques (an approach that was already 
considered in [5]), Fast Multipole Method, etc.; issues associated with such acceleration, however, 
will not be considered in the present paper. 

Our theoretical treatment relies on use of both existing and new analytical tools. In a first 
key step of the analysis the problem is re-expressed as a system of integral equations in a space of 
periodic functions. This is accomplished by means of yet another set of local cut-off functions, whose 
presence does not affect the actual numerical scheme. The use of periodic Sobolev spaces allows 
us to take full advantage of numerous results for approximation of Fourier series and interpolation 
operators [19j . as well as the theory of collective compact operators by Anselone [2]. Recasting the 
numerical scheme of quadrature type as a discretization method in L^-type (Sobolev) spaces gives 
rise to a number of difficulties. In particular, for the sake of the analysis we introduce Fourier 
projection operators, bi-periodic trigonometric interpolation operators, and a discrete operator [9j 
that produces a linear combination of Dirac delta distributions on a uniform grid from a smooth 
function input. The convergence analysis for the operators arising from the regular part of the 
original boundary integral operator follows the lines of the theory [9l [12] on periodic integral 
equations in one variable. The final (rather technical) element of our convergence proof is a detailed 
analysis of the integration error arising from the numerical polar coordinate integration algorithm 
for products of smooth functions, "shrinking" floating cut-off functions and bi-variate trigonometric 
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polynomials, in terms of Sobolev norms of the latter polynomials. 

This paper is structured as follows. In Section [2] we describe, in a compact form, the Nystrom 
method under consideration, and we state the main convergence results of this paper. In Section 
[3] we then recast both the continuous and discrete problems as systems of equations in spaces of 
biperiodic functions, we derive bounds, on various norms, of the main continuous integral operators 
in our periodic formulation, we establish unique solvability of the continuous system of periodic 
integral equations as well as the equivalence of this system to the original BIE, and, finally, we 
state the main approximation results in the biperiodic frame: norm convergence of the discrete 
operators to the continuous ones together with corresponding error estimates. In Section |4] we 
present key estimates on the singular operators that appear in the biperiodic formulation, and in 
Section [SI in turn, we provide the proofs of the main results stated in Section 3 and Section [2] — in 
that order. In Section [6] we describe and analyze a slight variant of the numerical method, where 
one of the interpolation processes is substituted by polynomial interpolation. Appendix [Aj finally, 
is devoted to the Sobolev error analysis mentioned above for the polar-integration of products of 
smooth functions, cut-offs and trigonometric polynomials. 

We conclude our introduction with two remarks concerning notation. 

Remark 1.1 To make a clearly visible distinction between points on the surface S and coordinates 
for their parametrization, we use boldface letters (e.g. r) for points on the scattering surface S C M^, 
and underlined letters (such as u) for points in M?. The coordinates of such points will be denoted 
according to u= (^1,^2). 

Remark 1.2 Throughout this paper the letter C denotes a positive constant independent of the 
parameters h and 5 and any other variable quantities appearing in the equation. When necessary 
a subscripted letter C is used — either to avoid confusion or to explicitly signify dependence on 
parameters other than h and 5. 

2 The Nystrom method 

2.1 The Boundary Integral Equation 

We consider the problem of time-harmonic acoustic scattering by a sound-soft obstacle with smooth 
boundary S in three-dimensional space: 

AU + K^U = inext(5), 

U = -U'""^ on S, (2.1) 

drU — ikU = o{r~^) as r —7- 00. 

Here k > is the wave number, C/™'^ is the incident wave, r := |x| and dr denotes the radial 
derivative. Letting denote the fundamental solution of the Helmholtz equation, 

a>.(r,r') ,^«-P(-|r-r'l) 



then the solution U of (12. ID can be expressed as the combined (or Brakhage- Werner [3j) potential 



U{r) 



du{r 



!;^^(r')dr' - i?] [ ^Jy,y'U{y')At\ 
) Js 



3 



where denotes the outward normal derivative on S and > is a coupling parameter. The 
density ■0 is the unique solution of the integral equation 

i^(r) + Aip{r) = -U'^'^ir) Vr € S, (2.2) 

where 

A^JJ{T) = ^i^(r,r')^(r')dr' := ^ {^^^^^ " ir/cl>«(r, r')) V(r')dr'. (2.3) 

Various choices of the coupling parameter rj have been proposed for accuracy and numerical effi- 
ciency; see e.g. [TJ [lOl [TTJ [5]. Note that the kernel K can be expressed in the form K = Kq + Ki, 
where 

/ sin(«;|r - r'l) ./ , sin(K|r - r'|) \ (r - r') • ly(r') 

^0 r,r' :=??^V ^-lUcos Kr-r' ^ ' , ,2 ^.4 

47r r — r V r — r / 47r r — r'r 



and 



, cos(K|r — r'l) K|r — r'l sin(K|r — r'l) + cos(K|r — r'l) (r — r') • i/'(r') 

i^i(r,r ) := -iry^^^— — ^ + 47r|r - r'| 47r|r - r'P ' ^^'^^ 



Remark 2.1 is easy to check that Kq E C°°{S x 5). T/ie kernel Ki, on the other hand, while 
not smooth, can he integrated with super- algebraic accuracy by means of a combination of a polar 
change of variables and the trapezoidal rule; see and equation (|2.23p . 

In this paper we focus on the Brakhage- Werner formulation presented above; a similar analysis 
can be used to treat the closely related Burton-Miller formulation 

2.2 Geometry 

The numerical method studied in this work relies heavily on the use of a system of local charts for 
description of the surface S. We will thus use a set of a number J of open overlapping coordinate 
patches {S^j^i that cover 5, 

J 

each one of which is the image of a smooth parametrization 

r^' : 

u:={ui,U2) I — > r^ui,U2), 

where C I2 ■= (0,1) x (0,1). We assume that can be extended to a smooth bijective 
diff'eomorphism between and so that, in particular, the Jacobians 

a^{u) := 

{a^ : D-' — )• M) are smooth functions of u. 



dr^{ui,U2) dr^{ui,U2) 

X ::: ^ U, 



dui 



du2 



J = l, 



(2.7) 



4 



The method requires exphcit use of a smooth partition of unity {uj^}j^i C C°°{S) on 5, subor- 
dinated to the covering (j2.6p . that is 

J 

oj-^ > 0, supp w-' C S-^ and ^^'^'^ = 1- 

i=i 

The hypotheses on availabihty of local charts and a partition of unity is not restrictive in practice: 
such parameterizations can be constructed for arbitrary geometries (see e.g. We will also 

assume that the boundary of supp a;-' is the finite union of Lipschitz arcs, a restriction which, 
again, is easy to accommodate [Ij. 
For any 6 > and r S 5 we let 

B{r,6) := {r' e S : |r' - r| < 6} 

and, selecting parameters < eo < ei < 1 that will otherwise be fixed throughout this paper, we 

define 

Si:= U B{r,2ei6)= [j B{r,2ei6). (2.8) 

r£suppa;J rGsuppoj-' 

Clearly there exists (Jq > such that for all j 



S(r,ei5o)nS^^^ =0 Vr G 5 \ 5^ (2.9) 

— as it can be checked by considering a pair of points that realize the distance between the bound- 
aries of supp a;-' and . In particular, this implies that 5*2^^^ C . The final element in our 
geometric constructions is a smooth function t; : 5 x [0, oo) — [0, 1] such that 

v{r, ■) = 1 in [0,eo(5o] 

Vr G S. 

v{r, • ) = in [ei5o,oo) 
Given 5 £ (0, Jq] we now define the functions tjs S x S ^ [0,1] 

%(r,r'):=^;(r,^^^^). (2.10) 

Clearly 



suppr/5(r, •) C S(r,ei(^), Vr € S. (2.11) 



In view of the definition of the sets Sg and (j2.9p we also have 

suppr/5(r, On 5^5/2 = VrG5\S^^. (2.12) 

Given 5 and considering the decomposition K = Kq + Ki in term of the kernel functions given 
in equations ()2.4p and ()2.5p , we define the regular part of the kernel of the integral operator (j2.3p 
as 

ir,eg,5(r,r') :=i^o(r,r') + (l-r?5(r,r'))/^i(r,r'). (2.13) 
Clearly i^rcg,(5 S C°^{S x S). The remainder is the singular part of the kernel, 

i^sing,5 ■■=K - K,eg,s = mKi, (2.14) 
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which, hke the kernel Ki, can be integrated accurately by means of a polar change of variables; see 
Remark 12.11 

We next introduce 

K'^{u,v) := K{r\u),r^{v))a^{v), 

(see equation ()2.7p ) with corresponding definitions of K^j^^^ and K^i^s (^^- equations ()2.13p 
and ()2.14p ). Noting that K*-', K^j^^^ and K^igg, are defined in x , we extend these functions 
by zero (possibly thereby introducing discontinuities on the boundary of D^) to the full product of 
squares I2 x /2- Clearly, 




Therefore, if ■0 is the exact solution of (|2.2p . then -0* := -0 o : — )• C {i = 1, . . . , J), is a solution 
of the system 

\/u£D\ i = l,...,J. (2.15) 



Remark 2.2 iVo<e i/iai the functions ip^ that appear in the integrals over I2 in equation (I2.15P 
are only defined in . However, they are multiplied by the cutoff function uj^ o which vanishes 
outside and thus provides a natural extension for the product throughout I2 ■ 

The solution of this system is used in Section [3] to reconstruct the solution of the original prob- 
lem (HJI). 

We can clearly distinguish two different types of integral operators in (12.15p . namely, integral 
operators with smooth and singular kernels. The discretizations of these operators are produced, 
accordingly, by means of two different strategies — as discussed in the following section. 

2.3 Discretization of the integral operators in equation f l2.15p 

For each patch , we select a positive integer A'^-, we let hj := l/Nj, and we introduce the grid 
points 

:= hjm = hj (mi, m2), < mi, 7712 < Nj — 1. 
We assume that these grids are quasiuniform: letting 

:= min A',-, h := 1/A^ = max hj, 
j=i,...,J j=i,--;J 

we assume there exists c > such that 

h<chj, j = l,...,J. (2.16) 
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Figure 1: Sketch of the the geometric layout of a single chart and its associated cut-off function. 

This is not a restrictive assumption in view of the assumed smoothness of S and C/™'^, and, therefore, 
of the solution (^■'). 

The algorithm under consideration is a Nystrom method which produces point-wise values of 
the unknown (tp^) at the discretization points xjn- The quadrature rules that ultimately define 
the method are described in what follows; for notational simplicity our approximate quadrature 
formulae use the set 

Ztv, := {m := (mi,m2) e : < mi,m2 < Nj - l} (2.17) 

of two-dimensional summation indices m. 

The approximate integration method used by the algorithm to treat the regular portion of the 
kernel is, simply, based on the trapezoidal rule, 

Ki^,si-^mv)dv^h] (2-18) 

— 3 

with the convention that K^lg s(' = DK Notice that for a smooth function ^ compactly 

supported in 

■= (rJ)-i(suppwJ) C D^, (2.19) 

(such as ^ = (u^ip) o r-' with ip G C°°{S)) the rule embodied in equation (|2.18p gives rise to 
super-algebraic convergence. 

To approximate integrals of the form 

/ i^sing,5(M,2£)e(l£)d^, eeC°°(/2), supple f]^ (2.20) 

that include the singular part of the kernel, the algorithm uses polar coordinates around the singu- 
larity. To properly account for contributions arising from various regions delineated by local charts 
and overlaps we let 

rji ■= (ri)-i o r' : (r^)-i(5^ n S^) C ^ 

s'J := n Si cs^n s^, -.= {vTHs^' n sj) = {vTHsf); ^^'^^^ 

see Figure [2] and equation (j2.8p . A close inspection of the definition of the sets Sjj shows that 
B(r, ei6) n suppw-' = if r ^5/2' Since, additionally, 

suppi^sing,5(n, •) C suppr?5(r^(n),r^(-)) C {r^y' {B{r^{u),ei6) D S^) , 



/ 

JI2 
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we conclude that the integral ()2.20|) with supp^ C , vanishes for all u such that r*(u) 

i.e., for u ^^/2- "^^^ algorithm therefore only produces approximations of ()2.20p for u £ ^^5/2- 
To do this consider the relations 

1 r27r / poo \ 

= - i^J Jp\Kl,,si^^L''{u) + Pm) ^ir^'{ll) + Pm)<ipjde (2.22) 

where e{6) := (cos 0, sin 0) and note that the function 

(P> G) ^ \P\K1,,S i'Xm^) + Pe{e)) (2.23) 

is smooth (as shown in Section HI cf. O |6l US]), 27r-periodic in 6 and compactly supported as a 
function of the variable p. In what follows we temporarily let z = (2:1,^2) := L^^iHi): so that the 
dependences on the patch indices (i, j) and parametric coordinate u are not displayed. 

The integral in the angular variable is approximated using a trapezoidal rule on a uniform 
partition of [0, 27r] in Qj subintervals of length kj = 27r/0j. Therefore, we consider the angles 

el:=pk„ p = o,...,e,-i 
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and approximate (j2.22p by 

y E /^^ig,5(M,l + />e(0^))IHCU + /^e(e^))dpj. (2.24) 

Recall that a uniform Cartesian grid in I2 with mesh length hj has been introduced in the approx- 
imation of the regular part (I2.18p . We will now use points of this grid to approximate the radial 
integrals in (j2.24p . 

If I cos el\ > V2/2, that is, modulo 27r, Ol G [-7r/4, 7r/4] U [37r/4, 57r/4], we look for the intersec- 
tions 

{x = z + pe{ei) : p eR}n{x = {qhj,0 : ^ G M}, q = 0, . . . , Nj - 1, 

i.e., the intersections of the double rays stemming from z with angle ±9p with the vertical lines of 
the uniform grid: the points of intersections are given by 

{qhj,Z2 + {qhj - zi) tan ei), q = 0,1, . . . , Nj - 1; (2.25) 

clearly, the distance between two consecutive points is hj/\ cos 9p\ < ^/2hj. Alternatively we can 
express these points in the form 

z + Pg'iz,9l)e{ei), where p,^(z,0) := i^^. 

For these values of 9p, the corresponding integral in (j2.24p is approximated by 

h,c{9l) ^ig,5(M,l + P^(£,^^)e(e^))l4^(l,^^)ie(l + P^(l,^^)e(e^)), (2.26) 

q=0 

where c{9) = 1/| cos 9\. For the remaining angles, intersections are computed with the corresponding 
horizontal lines: 

{x = z + pe{9l) : p£R}n{x = {^,qhj) : ^£R}, q = 0, . . . , Nj - I. (2.27) 

(The quadrature points are deliberately selected to lie on lines with one of the coordinates equal 
to qhj] the algorithmic reasons for such selections are discussed in Section [6l) Using the angle- 
dependent weight 

c(0) :=min{^-,— (2.28) 
I cos t'l I sm 9\ ) 

as well as the nodal points radii 

qh — zi 

p\{z,e) :=< 

I sin 6*1 

expression (j2.26p provides approximation of all the integrals in (j2.24p . 



-, if |cos6l| > \/2/2, 
I cos H\ 

-, otherwise, 
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In sum, we define our discrete singular operator (which depends on hj and kj, or, equivalently, 
on Nj and ©j) by 



p=0 q=0 



if n G il^-'^g ' 
otherwise, 



where ^^^^2 ~ ('^*) ^("^l^/a) ~ ('^*) ^("^^ ^h^)- what follows we use the scaling 



■'s/2 



(2.29) 



0.8 



0.6 



0.4 



0.2 
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Figure 3: Points used in the integration of the singular part. 



a > 0, 



(2.30) 



i.e., we assume that there exist C, c > such that for ah j, cA^- < Qj < CN-^". We point out 
that Pg{z_, ■ ) is discontinuous at {vr/4, 37r/4, 57r/4, 77r/4} but both side limits exist. It is immaterial 
which value (the right or left limit or even the average of both) is set as value to Pg{z, ■ ) at these 
points. 



2.4 The numerical method 



We are now in a position to lay down a fully discrete version of equations ()2.15p . The unknowns 
are taken to be approximate values 



(f^, ieQj^:= {rnGZ^. : r-'(x^) G int(suppcj-')}, j = l,...,J 



(2.31) 



of the true function values ^-'(a;^). In what follows we denote by (p^ G C{^}jJ the array whose 
i= (^1,^2) entry is for all indices We consider the set of indices 

^*N- ■= {m-= {mi,m2) G : -Nj/2 < mi,m2 < Nj/2}, (2.32) 
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(note that the cardinahty of TL*^. is #Z|^r^. = j^TLj^j = Nj), the space of trigonometric polynomials 

T]\f. := span | exp(27rim; ■ u) : m G ^at^. }, 
the interpolation operator : C{0,-j.^) — )■ T]\f. defined by the equations 



0, ££ In, \ 



and the vectors G with components ujj := {r^ (x''^)) . Using these notations, the discrete 

Nystrom equations for the system (j2.15p are given by 

J J 

yi + T.^' E Ki,A^i^2Q<^n, + T.Hh^^^^^^ (2.33) 

for G C(r2'^) and j = 1,. . . ,N, where the bynary operator denotes componentwise product of 
arrays: for example, (f^ is the array whose £-th entry is given by ujjip''^. Once these point 
values have been computed, the Nystrom method provides the reconstruction formula 

i,l = -2Y,(h] Ki,A-^^k)<< + ^t^M^' -2f/'"^°r\ (2.34) 

{i = 1, . . . , J). With these definitions, we clearly have 

Finally, the J parameter-space continuous functions : /2 — s- C can be assembled into a single 
continuous approximate solution V'/i(r) defined on S: 

iljh{r):= Y J{T)'4){{{r^)-^{r)), where X(r) := {j : r G suppo;^'}. (2.35) 

ieJ(r) 

The main convergence result of this paper can now be stated; its proof is given in Sections [3] 
through [5l The regularity estimates given in the statement of the theorem are expressed in terms 
of standard Sobolev norms on the surface S (see [1] or any standard text on Sobolev spaces). 

Remark 2.3 The notation 6 ^ means that 

ch^ <6< Ch^, 

for some constants c > and C > 0, independent of h and f3. Hereafter, the parameter /3 G (0, 1) 
will be taken to be fixed. Dependence of constants on this parameter will not be shown explicitly. 

Theorem 2.1 Let ip be the solution of ()2.2p . and assume 5 h^ with 13 G (0,1). Then for 
sufficiently large values of N the Nystrom equations (|2.33p admit a unique solution. Letting ip^ 
denote the reconstructed function on S as defined by (|2.34p and (|2.35p . we have the error estimate 

Uh-i^Wms) < Ct(|log/i|t/i*(i-^) (2.36) 
for all t > 2 (see Remark \l.S\) . Finally, for all e > and t > 2 

||V'/.-^||Lo.(s)<Ci,,(/i*(i-^)-(i+^)/^ + /i*-i)||V^||H*(5)- (2.37) 



11 



Theorem 12.11 tells us that for a smooth surface S and a smooth right-hand side U^^^ (from which it 
follows that ^ S C°°{S)) the Nystrom algorithm under consideration converges super-algebraically 
fast. 

3 Biperiodic framework 
3.1 Continuous equations 

The analysis of the method will be carried out by recasting the system (j2.15p as a system of 
periodic integral equations with all unknown functions defined on the unit square l2- To introduce 
our periodic formulation we utilize a second family of cut-off functions, a3^ G C°°{l2), that depends 
on 6 and satisfies the following assumptions: 

suppC^ C = (r^rHSl,,^), (3.1a) 



35/2 V ^ ^'^3(5/2 



1 in = {r^)-\Sl) D , (3.1b) 
\\dj:>l\\oo,h < Ca^"'-' Va > 0, (3.1c) 
where for a given non-negative bi-index a = {ai, 02), 

dd{y^) = dZidZii{ui,u2) 

and II • 1 1 00, 72 is the L°°{l2) norm. We will use the characteristic function 

'"o(^) ■■= [ otherwise, ^^'^^ 

which can be viewed as the limit of S;^ as 5 — )• 0. 

Using these functions, we define the following periodic integral operators 

A^ig/ := (^'or^)f K;^g,,(-,2i)5](2i)e(n)dn, (3.3) 
J I2 

Alg/ := (^^or^) /■ K^^g,(.,n)i}KM)?(^A)d^, (3.4) 

reg,<5 ' sing,<5 



h 

:= (Al, + Al .)e, (3.5) 



as well as the right-hand sides W := — (a;*C/™^) o r*, properly extended by zero to the full unit 
square l2- If u G /2 \ ^25' then r*(u) S^^ and 

supp(i^*/j^g,5(u, ■)oJs) C supp?75(r*(u), • ) n suppw^ C supp%(r*(u), • ) n fi^^^g = ^ 

by ()2.12p and ()3.1ap . which implies that 

(Alg/)(n)=0 yu€l2\ng. (3.6) 

Letting tp^ ;= ip o be as in equation ()2.15p and since uig [uj^ o y^) = oj^ o (see (I3.1bp and 
note that supp(a;-' or-') C ^^^"'), we see that the functions 

(t>i ■={u:^ij)or^ = {uj^ or^)^i : I2^C j = l,...,J (3.7) 
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extended by zero outside $7-' = supp(a;'' o r-') constitute a solution of the system 

\'i>' + E A:ig,.'^^' + E ^,s<t>^ = u\ z = 1, . . . , J. (3.8) 

The system (|3.8p can be considered as a system of equations for the vector {(j)^) G (L^(/2))^ for 

a given any right-hand side iV^) G (L^(/2))^- With this understanding, Theorem 13.61 shows that 

the system (j3.8p has a unique solution for any right-hand side (C/*) G (L^(/2))'^- (It follows that 
for the particular right-hand side (|3.8p the solution is (j3.7p — which, clearly, is independent of 5 in 
spite of the (^-dependence of the system of equations (j3.8p .) Theorem 13.51 shows that for the case 
:= —{oj'^U™'^) o r*, the solution of (|2.2p can be reconstructed from the solution of p.Sp . 



3.2 Discrete system 

In order to obtain a continuous system of equations from the fully discrete system (j2.33p we intro- 
duce an interpolation operator Qjv^. : C^{l2) '^Nj given by 

(Q7V,6(4) = ^(4) V^GZ^v,. (3.9) 
The discrete counterparts of the functions (j)^ are 

cPi:=ClM,{{u:^or^)i;{), (3.10) 

where ipj^ is defined by (I2.34p with (pjn obtained as the solution of the system (12.33p . Note that 

{0J-' LD-' if £ G fl-' 

0, otherwise. 

On the other hand, (w^(/>;^) (x;^) = ^oi all i, since (pj^ix-'g) = 0, for all i ^ and inside 

the functions do not have any influence by (I3.1bp . Therefore ()2.33|) can be equivalently 

re-expressed as an equation for continuous biperiodic functions {cl)\, . . . ,<p'l) such that for all i G 

{1,...,J} 



J 



j=l meZjVj 



14 + ^Q^,(/.,2 J2 (^^°r^)^;^g,5(-,3lk)^Ki£kK(^^ 

J 

+ Y.qN^{{u^or%l^4) = Q^^U\ (3.12) 



3=1 

Note that if the system (j3.12p is uniquely solvable, the solution belongs necessarily to TjVi x • • • x Tjv^ . 
The solutions of (|2.33p and (|3.12p are related by the formula (|3.1ip . 

For the sake of our analysis we recast the system (j3.12p in an equivalent operator form. To do 
this, we introduce the discrete operators 

Al, . . := (^^ o r^)Lt, (3-13) 
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and, for ^ G C°(/2), we define (cf. H and [T^) 

^N,i:=h] (3-14) 

— J 

where 5u denotes here the Dirac delta distribution at the point u G l2'- 6utp ■= '>P{u)- Further, 
we note that since the operators A*-^ ^ have continuous kernels, they may be applied to delta 
distributions: 

In view of these definitions, we clearly have 

— J 

and (I3.12P is equivalent to the equation 
J 

+ E {^,,s^^A + =Qn,U\ i = 1, . . . , J, (3.15) 

for the unknowns {(t)\, . . . , (p'l) G T^-^ x • • • x T^j. 

In order to recast this system of operator equations as a system defined in spaces, we insert 
the orthogonal projections Pat^. : L'^{l2) — ^ TTjv^. and write (I3.12p in the form 

J 

l^l + E (A:ig,.Div,P7V,<^i + A^.,,s,,^N,K} =Qn.U\ i = 1, . . . , J. (3.16) 

i=i 

for the unknowns {(p\, . . . , (/)j[) G L'^{h) x • • • x L'^ih)- 

Remark 3.1 iVote that the operator F]\f. of projection from to the the space of trigonometric 
polynomials make it possible to recast the fully discrete equation (j3.12|) (whose unknowns are ap- 
proximate point values of the continuous solution of equation (j2.15p ) via its version (j3.15p . posed 
in the space of trigonometric polynomials, as an equation (j3.16p in the complete space L^. 

3.3 Mapping properties of the operators introduced in Section 13.11 

As a result of the work in Sections 13.11 and 13. 2^ our overall integral equation and its Nystrom 
discretization have been re-expressed as equations (|3.8p and (j3.16p . respectively, which are posed in 
terms of unknowns defined (and compactly supported) on the open unit square I2. Such functions 
can naturally be extended to period- 1 biperiodic functions defined in M?, on which our analysis is 
based. In this section, in particular, we study the mapping properties of the operators introduced 
in Section [3.H when viewed as operators defined on spaces of periodic functions. 

To do this we first consider the biperiodic Sobolev spaces, see e.g. [19j, and we study their 
connections (that result through our use of local charts and partitions of unity) with the classical 
Sobolev spaces on the surface S. The Fourier coefficients of any locally integrable biperiodic function 
^ : — )• C are given by 

^{m) := / S,{u)emiu)'iu, em,(2i) 6xp(27rim • u), m&7?. 
J12 ~ ~ 
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For arbitrary s G M, the Sobolev norm 

iieii. := (le(o)i^ + Ed^ii' + \^2\'min)\'y^' (3.17) 

rn^O 

is well defined for all in the space T := spanjcm : m £ Z^} = UtvTat of trigonometric polynomials. 
The Sobolev space is defined as the completion of T under the norm || • \\s. Note that is 
the space of biperiodic extensions of functions in L^(/2). By a simple density argument, we can 
define the Fourier coefficients for any element of and any s < 0. It is also possible to define 
partial derivatives of any order: for a given non-negative bi-index a = (ai, 02), using the notation 
\a\ = ai + Q!2, we see that da is a bounded operator from to I for all s. 

The atlas introduced in Section [2] for representation of the surface S together with the H'^ norm 
just defined gives rise to a definition of the Sobolev norm 




(3.18) 



for any U E C°°{S) and any s G M. (In this formula it has been implicitly assumed that, even 
though r* is only defined on D'' C I2, the function (u^U) o r* can be extended by zero to a function 
in C°°(/2) — since the support of o r* is contained in — and the result can be extended as a 
smooth biperiodic function to all of M?.) The space H'^{S) can then be then defined, for instance, as 
the completion of C°°{S) in the above norm: this definition is equivalent to the classical definitions 
given in e.g. [I], [H §1-3.3] and [18, Chapter 3]. 

Lemma 3.1 For all s and j S {1, . . . , J} 

\\{uioo{r^r'\\Hs(s)<Cs6-\'\ms veGi^^ vj, (3.19) 

where Cs is a constant independent of 6. 

Proof. Because of ()3.18p we only need to prove that the maps 

e ^ Tf^ := {oj' o r'){^i o r^*) (^ o r^*) (3.20) 

V ' 

= Ss 

map -ff* to for every s and that their norms are bounded by a multiple of (5^I*L To do this we 
first note that 

svipp0^;c{rTHS'nS'), (3.21) 
(which is the domain where r-'* is defined), and that, in view of equation (j3.1cp . we have 

\\dag^\\oo,I,<Cad~^''\. (3.22) 

Since the operator T;J* is given by multiplication by the C°° function ffg' preceded by application 
of the smooth (J-independent) diffeomorphism r-'* : {r^)~^{S^ Pi 5*) — )• (r'')~^(5-' H 5"*), using the 
differential form of the Sobolev norms for positive integers s together with equations (j3.2ip 
and (|3.22p . we see that 

\\Ti'^\\s<Cs5~^ms yUH'^ s = 0,1,2,... 
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Together with the interpolation properties [19] of the Sobolev spaces i/^, this bound estabhshes 
the result for all s > 0. The result for s < follows from a duality argument. □ 

Here and in the sequel, || A||x~!.y denotes the operator norm of the bounded operator A : X ^ Y 
between the Hilbert spaces X and Y. 

Proposition 3.2 For all s ^ and all indices i,j, the integral operators A^-' : if* — t- H'^^^ given 
by equation (j3.5p are continuous and we have 



Proof. We can write A^'' as the composition (from right to left) of the maps i— t- (a;^^)o(r-' )^^ with 
A (the combined integral operator given in ()2.3p ) and then (p i— )• (wV) ° ^ - Lemma l3.ll provides 
a bound for the norm of first of these maps, whereas the norm of the last function as a map from 
H^{S) to the periodic Sobolev space is clearly bounded in view of the definition (j3.18p of the 
surface Sobolev norms. The result therefore follows from the fact that A is a bounded operator from 
H^{S) to H^'^^{S) for all s £ M — as it follows from standard results concerning pseudodifferential 
operators on smooth surfaces (see e.g. Chapters 6-9]). □ 

The next result studies the mapping properties of the regular and singular part of A^-' in the 
frame of the periodic Sobolev spaces and gives some estimates for the continuity constants in terms 
of S. In some of the arguments, it is convenient to use the space T> := DseRH^, consisting of smooth 
biperiodic functions. 

Proposition 3.3 For all i, j = 1, . . . , J , and s, t £ M the operators 

are continuous. Moreover, for all t > > s and all r £M there exist positive constants Cs,t <ind Cr 
such that 

WA'i^jH^-.m < Cs,t6'~\ (3.23) 
\\^Lj\Hr^Hr^^ < i-^'^+^'i>. (3.24) 

Finally 

W^IsJho^ho < Co. (3.25) 

Proof. Let 

F:i^^si^,v) ■■= {J o v%u)Kl^^,{u,vp,{v) 

be the integral kernel of the operator A^^^ ^. This function is well defined on x and can be 
extended by zero to the rest of I2 x I2 thanks to the cut-off functions that appear in its definition. In 
particular, Fli^ r admits a smooth biperiodic extension to x M^, and, thus A';-'^ is a continuous 

operator from to for all s, t G R (see [19^ Theorem 6.1.1]). Since A^j^^g ^ = A^-' — A*-^g ^, using 

Proposition 13.21 the claimed mapping properties of A*-?^g ^ follow directly. 
To establish p.23p we consider the operators 

■■= + := D^^i + ^(0). (3.26) 
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(Note that : — )■ is a bounded isomorphism for all s.) Letting 

G{n,v) := D^F^^,{u,v) + [ Fi'^^,{w,v) dw, 

we note that G £ C°°{M? x M?), that G is periodic in each of its variables, and that, when restricted 
to /2 X G has compact support. Now, for all r/ G 2? we have 

AmA:^eg,5Anr/ = / G{;v)Dl^{v)dv + m f G(-,v)dv 

D'^G{-,v}7]{v)dv + 9i{0) [ G{-,v}dv. 
h J12 

Further, bounds of the form 

sup \D^DlFii^^,{u,v)\ < C^,„r2™-2", (3.27) 

follow from the definition the kernel functions -f^^eg 5 (^^^ (|2.13p and subsequent lines) and from the 
assumptions (IS.lcP on S;^. It follows that, for ^ E 2? we have 

■2m— 2n II ^ ii 

-2n, 



where we have used the integro-differential form of A^A^-^g ^A„ and (j3.27p . This inequality estab- 
lishes ()3.23p in the particular case t = 2m, s = —2n; the result for general t > > s then follows 
by interpolation [19j. 

To establish (|3.24p . in turn, we first note that (I3.23P implies 



f \\Al^jHr^HO<Cr,oS', ifr<-l, 

\\A%^jHr-,Hr+^ < Cr,r+ld-\ if -1 < r < 0, 

||Aj^ .||^o_H^+i < Co,.+i<5-('^+i), if r > 0, 



or, equivalent ly 



\ \ij II ^ ^/ s:- max{-r',r+l,l} 



Since A^j^^g g = A^^ — A*-^g g, and in view of Proposition 13.21 we have 

WKinsM^Hr + l < WA'iWHr-.Hr + l + 1 1 A:^^^, 1 1 ^^. + 1 < C,rH + C'^d' 

< max{-r,r+l,l} 

which establishes (I3.24p . 

Inequality (j3.25p , finally, follows from (|3.23p with s = t = and Proposition 13. 2t 

WA'sWho^HO < ||A'/||^^o_^j^i < C, 

with C independent of S. □ 

The limiting case 5 = is studied next. 



17 



Proposition 3.4 For all the operators Aq : — )• are bounded. 

Proof. Clearly, Aq ^ = A^-J^ (^^oO) where uJq coincides with the characteristic function of set fi-' (see 
(j3.2p ). Since the operator of multiplication by Wq defines a continuous operator in H^, the result 
follows from Proposition 13.21 □ 



3.4 Analysis of the continuous system 

Theorem 3.5 Let (cj)^, . . . ,4"^) be a solution of (|3.8p with [/* := —{uj^U^^'^) or*. Then the solution 
of (|2.2p can be expressed in the form 

V^(r) := '^'■((r^rHr)), 

see equation (|2.35p . 

Proof. Because of the particular form of the right-hand side as well as the presence of the 
factor u)^ o r* in the operator A^-' (see equations (|3.3p . (j3.4p and (|3.8p ). it follows that suppi?^* C 
supp(a;* o r*) = Q,\ Therefore, by ()3.1bp . cl;^0* = 0* for all i. 
Consider now the functions ip^ : ^ C given by 




These functions are constructed so that 0* = (w* o r*)-;/;* for all i. Although the functions ip^' 
are smooth only up to the boundary of (the domain of the chart r*), the function can be 
smoothly extended by zero to M^. Similarly, we can consider the functions : S — )• C such that 
$*(r) = CiJ*(r)^/'*((r*)~^(r)) = (/)*((r*)~^(r)) if r G and are zero otherwise. These are smooth 
functions on the surface S and ip = X]/=i - 

Then, if r G S and i £ 2^(r), we can write r = r*(n*) with := (r*)~^(r), and 

J 

{l^P + Ai;){r) = 1 <A*(n*) + ^a;*(r)(AV')(r) 

jgX{r) j=l 




iex(r) i=i ieX(r) 
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Note that we have used the fact that ujgcp^ = 4P . This finishes the proof. □ 
The product space 



J times 

win be endowed with the product norm, also denoted by || • \\s- We can then consider the matrices 
of operators 

as well as the identity operator I : T-L^ — t- T-L^. With this notation, ()3.8p can be written in operator 
form as 

Bsct):={^l + As)(l) = \J (3.28) 

where 4> := {(f)^, cjP' , . . . , (f)"^) and 

U := ([/\ . . . , C/-^) = -(w-'C/'"") o r^'. (3.29) 



Remark 3.2 Note that, while the operator As depends on 5, equations (|3.3p - (|3.5p show that, for 
elements 4> = {(1)^, ■ ■ ■ ■,4''^) such that supp(/>* C supp(ci;*or*), As(t> is independent of 5. The following 
theorem shows that the operator in equation (j3.28p is invertible, and, thus, in view of Theorem \3.5l 
for right-hand sides of the form (|3.29p . the solution of equation (j3.28p is independent of 5 as well. 

Theorem 3.6 For all < 6 < Sq, the operators Bs = \Z + As : Ti^ ^ Ti^ are invertible. Moreover 

Proof. We will first show that Bs ■ T-L^ — T-L^ is invertible. Propositions 13.21 and 13.41 and the 
compact injection C prove that As '■ — l-L^ is compact (and uniformly bounded in 5). 
Therefore, Bs ■ "hP — )• T-L^ is bounded and Fredholm of index zero. Thus, it suffices to prove the 
injectivity of the operator. 

Let then 4> £ Ker^S^, that is, 

J 

i</.^ + J^Ay<A^=0, i = l,...,J. 

Arguing as in the proof of Theorem [331 it is clear that cD^f/'-' = for all j, which means that Ker 
is independent of 5. Since : — )• H^~^^ are continuous for all s, then cff G C°°{l2). Following 
the argumentation in the proof of Theorem 13.51 we define 




(3.30) 



so that = (w* o r^)ip'^ . We now set 

^(r):= ct>^{{rn-\r))=Y,^^{r), 
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where the functions : S ^ C are defined by extending ip^ o (r-' )"^ = c/)^ o (r-' )"^ by zero to 
S\ . Proceeding as before, we prove that ^■0 + = 'ind therefore tp = 0. 
From (j3.3Up we deduce that for uGD^, 

= -2V / K'^u,v) (l>> (v) dv = -2y2 [ K{r'{u),r')<^^r')dr' 

= -2 [ K{r\u),r')i;{r')dr' = 
Js 

since tp itself is zero and the functions are each supported in . Therefore (p^ = (w* o r^)^p^ = 
and the injectivity of Bs is proven. 

To prove the uniform boundedness of B^^ : %^ — ?■ %^ for > we proceed as in the proof of 
Theorem 10.9 in [16j. Recall that As ■ 'hP — s- T-L^ is uniformly bounded in 5 by Proposition 13.21 with 
s = 0. Since the injection TH} C 'hP is compact, the set {^5)5 turns out to be collectively compact. 
Moreover, it is easy to verify that lims-^o A^^, = .4o^ in for all ^ G Ti.^ . Applying that pointwise 
convergence is uniform on compact sets and the collective compactness of the set of operators {As} 
it follows (see [Ml Corollaries 10.5 and 10.8]) that 

11^0"^ A - A)A||wo-,wo ^ 0. (3.31) 

Consider now Cs ■= 2(X — Bq^As), which is uniformly bounded. Straightforward computations 
show that 

CsBs =Z + 2Bo\BoA5 - BgAs) =1 + 2BoH-^o - A5)A5. 

Therefore (j3.3ip proves that for 6 small enough CsBs is invertible with uniformly bounded inverse 
and therefore, since Cs is uniformly bounded, so is B7^. □ 



3.5 Convergence estimates in the biperiodic framework 

We now state a convergence theorem (Theorem 13. 8p in the biperiodic framework; as shown in 
Section [H this theorem is equivalent to our main convergence result. Theorem 12.11 The proofs of 
the two results presented in this section (Theorems 13.71 and 13. Sp . which require certain analytical 
tools that are developed in the next section, are given in section [5^21 
Let us introduce the operators 

Qh^, ■= (QNi^,!, ■ ■ ■ ,QnjCj) , 
Vh^ ■■= (PiVi6,---,P7V,0), 

Vhi := (D;vi6,---,D^^O), 

see Section [3^ and let 

Clearly, equation ()3.16p (and, thus, in view of Remark 13. H the main Nystrom system presented in 
equation (I2.33P ). can be re-expressed in the form 

Bs,h(l>h ■■= {\Z + QhAf^-DhVh + QhAfj^^Vh) cPh = QhV (3.32) 

where (ph '■= {4>h^ 4>hJ • • • > 4'h)- (Note that, although not explicit in the notation, the unique solution 
4>h of equation (j3.32p does depend on the parameter 6.) 
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Theorem 3.7 Let 6 ^ with /3 E (0, 1). Then 

lim \\Bs^h - BsW-^o^-^o = 0. 

Thus, for all h small enough, Bs,h '■ 'H^ 'H^ is invertible, with h-uniformly bounded inverse. 

Hence, for h small enough the numerical scheme admits a unique solution which depends con- 
tinuously on the right-hand side. 

Theorem 3.8 Let 5 ^ with (3 £ (0, 1) and let U be given by (|3.29p . Let 4> and 4>h be the 
respective solutions of Bs(l> = U and Bs^h4^h = Q/iU. Then, there exist constants Cg^t > for all 
t > s > and t > 1 such that 

u-(t^h\\s<Cs,th'^'-('^-'\mt. 

Note that cj> does not depend on 6 (see Remark I3.2p , although is a 5 dependent quantity. 

4 Estimates for the singular part of kernel and operator 

Lemma 4.1 The functions 

iu,p,0)^\p\K^{r\u),r^{r^\u)+pe{e))) 
are C°° in their domain of definition. 

Proof. The kernel function Ki can be decomposed (see (j2.5p ) as 

Fi(r,rO , F2(r,rO ^ 

where Fi,F2 G C°°{S x 5). Let 

fiiu,P,0) ■■= F,{r\u),r^iri\u) + pe{e))) £ = 1,2, 

gi{u,P,0) := p-V{u)-r^{r^\u) + pe{d))\' 

92iu,P,0) := p-^{r\u) - ir^\u) + pe{0))^ ■ u{r^\u) + pe{e)). 

It is clear that fi and /2 are C°°. Noticing that the functions 

z ^ \r'{u)-r^{r^'{u) + z)f (4.1) 
z ^ (r\u) - r^ir'^u) + z)) ■ i^{r^\u) + z) (4.2) 

satisfy the hypotheses of Lemma 14.21 it follows that gi and g2 are C°°. It is also clear that gi is 
positive for p ^ 0. The Hessian matrix at z = of the function defined in (j4.ip is the matrix with 
elements 

29.„^ (u) ■ r-'' (n) l<n,,n2<2. 
Clearly this matrix is positive semidefinite, and since its determinant is 

A\d,,r^{u)xd,,r^iu)\^ = Aa^iuf, 
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(see (j2.7p ). it is positive definite. Using the equality (|4.3p . it follows that gi{u,O,0) ^ 0. The 
function in the statement is 

-1/2 „ , -3/2 „ 
9l Jl+9l 12 92 

and the previous arguments show that it is smooth. □ 



Lemma 4.2 Let F be a C°° function in a neighborhood of the origin in M?. If F{0) = and 
VF(0) = 0, then the function 

f{p,9) :=p-'F{pe{9)) 

is C°°. Moreover, 

f{0,e) = ^e{9)-BF{0)e{e), (4.3) 
where HF(0) is the Hessian matrix of F at the origin. 

Proof. The result follows from an application of the Taylor formula 

h{p, 6) = h{0, 9) + pdph{0, 9)+p^ [ d'^h{p u, 9) (1 - u)du. 

Jo 

to the function h{p, 9) := F{pe{9)). □ 



The forthcoming analysis relies heavily on use of families {xs ■ < S < 5o} of functions for 
which there exist constants cq > and Ca for a > such that for all (5 {0 < S < 6q) 

XS G C°°(m2), (4.4a) 

X5ip,0) = Xs{p,9 + 27T) y{p,9)eR^ (4.4b) 

supp X5 C (-Co 6, Co 6) X M, (4.4c) 

WdaXsh^iR^) < Ca6~\^\ Va > 0. (4.4d) 

Proposition 4.3 For all u G ^2So ~ {^^)~~^{S'^ H Sg^^) and 5 < 6o, let us define: 

XuAp,0) ■■= ^u;*(r^(li))|p|i^ig,5(2^>i:''(2^) + P^W)^Ki:''(li)+/'^(^)) (4-5) 
Then the sequence {xu.s} satisfies conditions (14. 4p with cq and Ca independent of u. 

Proof. We start by considering the simpler functions 

^„,5(p,e) := \p\ Kll^/u,r!'{u) + pe{9))- (4-6) 

Since ^^^^ C {r^y^{S^ n 5^), the mapping r^^ := (r^)-'^ o : C h ^ h is weh defined. 

Therefore, the functions 

riu,5 (z) := r?5 (r* {u) , {r>' (u) + z)), 
fu{z) := \z\K,{r\u),r^{ri\u) + z)) 

are well defined in := {z : r^''{u) + z G D^} (that is = — l'^u))- 
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As a simple consequence of (j2.9p . if r G ^2^^, then B{r, ei6) C . Applying (j2.1ip and the fact 
that r*(n) G S2g^j, it follows that 



suppris{r\u), •) C S(r^(n), ei<^) C 5^. (4.7) 

Consequently 



suppr?„,5 C {r^)-\Bir'iu,ei5)))-r^\u) C {r^)~HS^) - r^\u) = 

Therefore supp rj^^s is strictly contained in fi^ and r]u_s can be extended by zero to a function in 
C~(]R2). 

On the other hand, for z G supp rj^^s we have 

111 = \r^' iu) + z- r^' (u) \<Cj\ {r>' (n) + z) - {r^' (u) ) \ 
= Cj\r^{r^\u) + z)- r'{u)\ < CjCid 

(we have used ()4.7p in the final inequality). Therefore we can fix c > independent of the particular 
chart number j so that supprju^s C B{0,c6) for all u G ^25o' Moreover, from (j2.10p we see that 

||9a%,<5||Loo(R2) < Ca5~l-I Vu G 0,'^^^. 

Thus, 'i]u,5ipe{9)) satisfies the conditions (j4.4p . 

Finally, since ijju,sip,d) = { fuVu,s){p cos 9, p sin 6) and since, by Lemma HT| fu[pe{9)) is smooth 
on the support of T]u,siP&i^))^ follows that {ipu^s} satisfies ()4.4p with cq and independent of 
u. In view of the inequalities (13. 1|) the same is true of the family {xu,5}, and the result follows. □ 



Proposition 4.4 There exits C > suc/i t/iai /or all N 'N and all S, < 5 < Sq, 

\\^Ls,5^n\\oo,I, < C\ log7V|3/2||^^||o V^jv G Tn. 

Proof. By ()3.6p . it suffices to bound the values of (A*-?^^^ s^^)(^) ^'-'^ — ^ ^2S- "^^^ polar coordinate 
form (I2.22P gives, for u G Q^'' 



2<5' 



(A:L5^^)(^) = / / XuAP,(^)^Nir''iyd + Pemdpd9 



-coS JO 



/cqS r27T 
/ Xu,5{p^^)^rnipei9))dpd9, 



where x«.5 is given by (14. 5p . Therefore, by the Cauchy-Schwarz inequality we have 

IIAlg/Tvlloo,/. < ll^ivllof E ^L^)'^'' (4-^ 



where 



fCo5 
-co (5 



-Dm,5 := max / / X2i,5(/0, 6*) exp(i/9m • 6(6*)) dp 



d0. 
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As a direct consequence of Proposition 14.31 it follows that 



co<5 



Co 5 



Xn,<5(p, 0) exp(i/jm • e{9)) dp 



< C5, 



and, integrating by parts, that 



Co (5 



Xu,5{p, 0) exp(ipm • §{9)) dp 



. ... exp(ipm-e(g)) ^ 

Co 5 



m • e{6) 



< C- 



m ■ ei9) \ ' 



where the constant in both inequalities is independent of u. Therefore 



27r 

Drr,. f,<C i min 



\os.N 



^ I mi cos 6^ + m2 sin 6^1 J y^mf + m| 



VmGZjv\{0}, 



where we have applied Lemma 14.51 for the last bound. Inserting this bound in ()4.8p and using the 
fact that 



E 



, <ClogiV, 



(this can be proved by comparison with the integral of l/(x + y )) the result follows readily. □ 



Lemma 4.5 For all m^l?\ {0} 

r 

/ min < 1, 



1 



mi cos + m2 sin ( 



dO < 



27r(l + log y/m'i + 



ms 



Proof. With some simple trigonometric arguments we prove that 

r>27r />27r 



I min{l, |mi cos^ + m2 sin0|~^}d0 = / min{l, |Re((mi — im2) exp(i0))|^^}d0 
Jo Jo 

= 4 / min{l, (wm^ + mg sin0) }dO. 



If we now shorten c := rnf + m'2, we can estimate 

rTT/2 pTT/2 



/ min{l,(c sin6')-^}d6l = c"^ / min{c, (sin 6')-^}d6l 
Jo Jo 

_1 TT IT 

~ ^ Jo 2^-^^^^ 2^^^"^^°^*^^' 

which finishes the proof. 



□ 
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5 Proofs of the main results 

5.1 Auxiliary approximation properties 

The operator P^v : U^i?* — )• T^v that cuts off the tail of the Fourier series and at the same time 
gives the best approximation in T^r for all s is: 

It is straightforward to check that 

\^Ni-i\\s<{'ih)'~'\\i\\t y^eH\ s<t. (5.1) 

The interpolation operators introduced in (j3.9p satisfy (cf. \19\ Theorem 8.5.3]) 

\\QN^-C\\s<Cs,th'-'m\t y^€H\ 0<s<t, t>l. (5.2) 

The following lemma studies the uniform boundedness of Qn as an operator from the space of 
continuous bivariate periodic functions to H^. (See [17\ Lemma 11.5] for a different proof of this 
result in the univariate case.) 

Lemma 5.1 For all N 

1 / \^/^ 

IIQjVello = ]v( E l^(^n)r ) ^ ll^ll-.^^ ^ C°(^2). (5.3) 

Proof. Given £ "^N we can write 

CN{^n) = ^ $iv(ini)exp(^ (n • 221)), n G Zat, 

and therefore, using the 2-dimensional discrete finite Fourier transform, 

When = Qat^, this gives the equality in (j5.3p . whereas the inequality is a simple consequence 
of the fact that #Zn = N^. □ 

We finally give an approximation result for the operator Djv defined in ()3.14p : 
Lemma 5.2 There exists Cg^t > independent of h such that 

\\GN'PNC-as<Cs,th'~'\mt veeF*, s<t<o, s<-i. 

Proof. The result is proven in ^ Lemma 6] (see also [12J ) for univariate 1-periodic functions (with 
the restriction s < —1/2 instead). The proofs can be easily adapted to the current case. □ 
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5.2 Proofs of the results of Section 13.51 

Recall that N has been taken to be the largest of the discrete parameters Nj and therefore hj < h 
for all j = 1, . . . , J. Also, by (|2.16p . we can bound hj^ < Ch"^ and | log hj\ < C\ log h\ whenever 
needed. 



Proposition 5.3 For all t > 2, there exists Ct > such that 

Proof. Let us first consider the decomposition 

^Ls-^N.^i^^s^N,^N, = A:^g,(I-D^^,P;v,) + (I-Q7vJA;^g,(D^T^^,-I) 

+ (I-Q^JA:ig_,. 

By Proposition 13.31 and Lemma 15.21 for all t > 2, 

- D^^.P;v,)ello < Ct6-'\\{1 - D^,Piv,)CII-t < QVr*||C||o. 
The second term is bounded using ()5.2p . Proposition 13.31 and Lemma 15.21 

11(1 - QivjA:^g^,(D;v,pjv, - mio < cth'^'\\^i^,si^N,PN, - mt/2 

< Cih'/^5-'\\{BN,PN, - m-t/2 < C'/6-'h'mo. 
Finally, by ()5.2p . Proposition 13.31 and Lemma 15.21 

||(I-QivJA;^g/||o < a/.iA:^g/||, < c;/i*r*||C||o, 



and the proof is finished. 



□ 



Proposition 5.4 There exists C > such that for all ^at^ G Tw. and 6 > h 



sup 



sup 
sup 



i^,,s - ^,,s,k)^N,(M)\ < c[6-'h,\logh,\+5)u^^\\o, 



King,5-AsingAh)^^.(^) 



(5.4) 
(5.5) 
(5.6) 



Proof Recalling the definition of the discrete operator Agj^^g^^ = (w* o r*)Ly^ (the operator L 
is defined in ()2.29p ). it is clear that 



-'6,h 



5/2- 



(5.7) 



Using also (|3.6p . it is obvious that both operators in (j5.6p vanish for u£l2\ ^26- Also (j5.5p is a 
simple consequence of ()5.7p and Proposition 14.41 



26 



For u G r^^-'^2' recalling that ^Njip, 6) = ^.N^ir'^iu) + pe{0)), 

/OO f'27V ^ ^ 

/ {Xu,5CN,){p,0)dpde - Qh,k,f{Xu,5CN,), 
-OO Jo 

where Xu,5 is defined by (j4.5p and the quadrature rule Qh,kr( is given in ()A.2p . By Proposition 14.31 
we can apply Theorem lA.l^ which estimates the above quadrature error in terms of the constants 
that appear in ()4.4p . Since these constants can be taken to be independent of u E (this is part 
of the assertion of Proposition I4.3P . then (j5.4p follows readily. □ 

In the following sequence of results we will use the geometric cut-off operator 
= (GiCi, . . . , GjO) := ((1 - ^DCi, ■ ■ ■ , (1 - ^0 
where, as a reminder, uJq is the characteristic function of the domain . 
Proposition 5.5 For all t >2, there exists Ct > such that for all (^^^ S Tf^f. and 5 > h 

IKAlg,, - Al,,,,JC^,||oo,/. < aS'-'h'MNMt + C5\\G,CNMh- (5-8) 

Proof. Following the argument of the preceding proof, but using ()A.4p of Theorem lA.H it follows 
that for any integer t >2, 



sup 



(Al,. - Al,,,)e^,.(n) < a5'-'h'Ur,^\U. (5.9) 



From the definition (|2.8p . it follows that if r then B{r,ei6) n suppw-' = 0. Therefore, if 

snvpKii^/u, ■)nn^ C supp%(rHn), • )) n 

C (r^r^(S(r»(n),ei(5) nsuppw-') = 0, 

where we have used (|2.1ip . This means that if u G 1^25 \ ^5/2 ' o^ly the value of on /2 \ (where 
(,Nj = is relevant. We then change to polar coordinates as in the proof of Proposition 14.41 

to obtain: 

/cqS f'27r 
/ IXM,5(p,6')|d/)d6' 
-C06 Jo 

< C6\\GjCN,\\oo,h- (5.10) 

(In the last inequality we have applied Proposition 14.31 according to which the constants cq and C 
do not depend on u.) Equation ()5.8p now follows from equations ()5.6p . ()5.9p and (jS.lOp . □ 
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Proposition 5.6 For all e > there exists > independent of h and 6 > h such that 

\Ar' - ^4^,^^, < Ce {s-'-'h + 5-'h\ log h\ + log . (5.11) 

Furthermore, for any integer r > 2, 

II {Af^ - Q/.4fn)^||o < CriS-'-h^UWr + S\\Vh^ - ^lloo,/, + ^ 11^^ ||oo,/,) G (5.12) 

Proof. We start with the decomposition 



Al,.(I-P^) + (I-Qiv.)Al,.P 



(5.13) 



For all t > 1, Proposition 13.31 with r = — 1 and (|5.ip imply 

iiAlg,,(i - Piv,)eiio < - PN,a-i < (5.14) 

For i > 1 , by (j5.2p , Proposition 13.31 with r = t — 1 and the inverse inequality ()A.23P , we can bound 

IKl-Q^JAlg^^P^^-ello < Cth'\\A'i^^^,FNMt<C'tS-'h'\\FNMt-i (5-15) 

< C;r*/i||P^^CI|o < C[5-'hU\\o. (5.16) 

For fixed ^ we introduce the auxiliary function g := (A^j^^g ^ — A^j^^g g /i)Pa^j?) ^^^^ by Lemma [5?T| 



(5.17) 



We can then split I^Ni as follows: 

A key point is the fact that i^Bjy- < C5Nf, which is proved in Lemma l5.7i From (j5.6p it follows 
that g{x]^) = for every m G j4iv-. Also, by (|5.4p 

(/i? E Ig&^y^' <c[6-'h,\logh,\+6)\\FN,^\\o. 

meCNi 

Finally, from (j5.5p . we obtain 



0- 



Going back to (j5.17p . we have proved that 

IIQ^.(Alg,, - Aig,,,,)P^,eilo < c(6-'h\ log h\ + log hf') \\FnMo- (5-18) 
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A direct estimate from (|5.17p . using (jS.Sp now, gives 



IIQ^.(Alg,5- A:j,g^,^,)P;v,ello < a^'-'^/iHlellr + C5||G,P^^^|U,/,. (5.19) 

Using the decomposition (j5.13p and the bounds ()5.14p with t = 1, ()5.16p with t = 1 + e and 
(|5.18p we can easily prove (j5.1ip . At the same time, (|5.12p follows from the same decomposition, 
using now (j5.14p with t = r + 1, (jS.lSp with t = r and (j5.19p . In a last step we use the fact that 
the cut-off operators G-' are bounded in L°°(/2). This finishes the proof of the Proposition. □ 



Lemma 5.7 There exists C independent of of 6 > h such that 

#{m G Z^, : e n% \ J^^^J < C5NI 

Proof. We start by setting := d{suppu^), which has been assumed to be the finite union of 
Lipschitz arcs. Therefore, for every 6 > 0, we can pick up a finite set C F-', with the following 
properties 

F C U B{r,ei6) with #F^ < CS'^. 

The constant C that controls the size of the discrete set F;^ depends on the Lipschitz constants 
related to F-'. On the other hand, for 6 small enough 

5^,\(suppa;^)C \jB{r,4ei6) 
rer 

and therefore 

Sis \ C Sis \ (supp ^^■) C U B{r,5e,S). (5.20) 

We now go back to the parametric domain using (r*)~^ to define the sets 

Or := {r')-\S' n B{r,5ei6)) r € F^^. 

Since ju — v| < C^\r^{u) — r^{v)\ for all u,vG D^, each set Or is either empty or can be surrounded by 
a closed ball of radius 5C*ei5. This gives a collection of points {uf, . . . C with ks < C6~^ 

such that 

^s\^/2 = i^T\sg\s'^^,) = {rT\s^n{sis\sl/,)) 

ks 

C \JOrC[J{v£R^:\v-]4\<5C'ei5}. 

The number of points of the uniform grid {x^ : m G Zat-} that fit in a closed disk of radius r 
is bounded by {2rNi + 1)^. Therefore, the cardinal of the intersection of the uniform grid with 
^25 \ ^5/2 '^^^ bounded by 

ks{lQC'ei6N, + if < C5-^{l{)&5Ni + 1)^ < C'SNf, 

where we have applied that 1 < 5Ni (which is implied by 5 > /i), and the result follows. □ 
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Proof of Theorem 13.71 This result follows now from an adequate choice of parameters in 
Propositions 15.31 and 15.61 Since S « /i^, then 5-^h ^ h^'^. We choose e = {I - (3)/{2(3), so that 
6~^~^h ~ We then apply Proposition 15.31 with t = 3 and the first bound of Proposition 

15.61 with the above e. This yields the bound 

11^5,, - ^^ll^o^^o < + + h^-P\ log h\ + hPI^\ log 

which ensures convergence as /i — )• 0. The uniform bound for the inverses of Bs (Theorem 13. 6p 
ensures then a uniform bound for the inverses of the operators B^^h when h is small enough and 
5^hP. ' □ 



Proof of Theorem [Ml For ip £ H''{S), let us define 4> := ((^V) o r\ • • • , (^'^'0) ° r-^) e W. 
By construction, Q(f) = 0. Applying Proposition 15.31 and the second bound of Proposition 15.61 it 
follows that for any integer r > 3 

\\iBs,h - Bs)H0 < Crih'-^^'^^imr + h^WVhCt^ - Hoc,!,). (5.21) 

Taking now 

r > max{l//3,/3/(l - /3)}, (5.22) 

so that 1+/3 < H-r(l — /3) < r, from Sobolev's embedding theorem and the approximation estimate 
(|5.ip we obtain 

h^\\Vhcl>-cl>\\ooj,<C^h''\\Vhcl>-cj>\\,+^<Cr,^h^'^^^^^ 

Applying this estimate in (|5.2ip we deduce that for every integer r > 3 satisfying (|5.22p . we have 

\m,h-Bs)(t>h<Crh'-^'-''^\mr. (5.23) 

Since /? is fixed, the dependence of the various bounding constants on the parameter /3 is dropped 

from the notation in what follows (recall Remark 12. 3p . 

For ip £ H^{S) we now define M.s,hip '■= {^5,h — Bs)4> with as above; clearly Ms,h '■ H'^{S) — ?■ 
is a continuous map. In view of ()3.18p equation ()5.23p can be re-expressed in the form 

\\Ms,hHo < Crh^^'-^^miHr^s) E H'^iS). (5.24) 

Since \\Bs^h — -^^II-ho-^W" ~^ 0' it follows that the sequence Ais,h is uniformly bounded, 

\\Ms,hi^\\o < Cmno^s), 

and therefore, using Sobolev interpolation theory [181 Appendix B] for the operator M.s,h, it follows 
that (I523D, and therefore (fOS]) . hold for ah r > 0. 

Let now cj> and cj)h be the respective solutions of Bscj> = U and Bs^h<Ph = QftU. From Theo- 
rem [321 it follows that for h small enough, the inverse of Bg^h ■ T~i^ T~i^ is uniformly bounded with 
respect to h. Hence 

U-Mo < C\\BsM<t>-(t>h)\\o<C\\Bs<t)-BsMo + C\\{Bs,h-Bs)ct)\\o 

= C\\IJ - QhVWo + C\\{Bs,h - Bs)(t>h. (5.25) 

But, from ()5.2p . it follows that for t > 1 

||U - Q;,U||o < Cth'\\\J\\t = Cth'WBsMt < C'th'UWt. (5.26) 
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Thus, substituting (|5.23p and (|5.26p in (j5.25p . we conclude that for all t > 1 there exists Ct > 
such that 

U-(t^h\\o<Cth'('~''^U\\t. (5.27) 

Error estimates in stronger Sobolev norms can now be obtained by means of the inverse inequali- 
ties ()A.23p together with the fact that Vh provides the best approximation, for all periodic Sobolev 
norms, in the discrete subspaces of trigonometric polynomials. Indeed, for s > we have 

110 - Ms < 110 - Vh(t>\\s + Ch-'\\Vh^ - 0h||o < 110 - Vh^s + Ch-'U - MW (5.28) 
The proof now follows by substituting (jS.ip and ()5.27p in ()5.28p . □ 

5.3 Proof of Theorem 

In view of the reconstruction formula of Theorem 13.51 and since Bs4> = U, the exact solution of 
(j2.2p can be expressed in the form 

^(r) = (a;V) ((rr ^r)) = J] .^^ ((rr ^r)) 

ieX(r) iGX(r) 




Similarly, for the discrete Nystrom solution we have the relation 
V^,(r) = (a;Vl)((r^)''(r)) 

iGX(r) 

= E - E ((r^)-^(r)) , 

iGX(r) j=l i=l 

which follows by re-expressing equations (j2.34p and (j2.35p in terms of (/>;^ = Qat^. ((w-' o r^)7pj^). 
Since ||V'||_f/=(s) = Il0l|s) by Propositions 15.81 and 15.91 below we thus have 

< C,(h''^ + |log/i|3/2/^*(i-^))||0||, = C,[h'-^ + |log/i|^/^/.*(^-^))||Vlk*(5), 

and (|2.36p follows. Note that the constants Ct depend on j3 (see Remark l2.3p . The L°°{S) bound 
(j2.37p follows from similar arguments, and the proof of Theorem 12.11 is thus complete. □ 

Proposition 5.8 For all t > 1 there exists Ct > such that 

\\^1,,6^nA - ^l,,6^h < Cth^'-'-^^mt. (5.29) 
In addition, for all e > and t > 2 there exists C^^t > such that 

\\^e,,S^NA - ^e,,5^\\oo,I. < C,,th'^''^'>~^'+^'>(' UWt. (5.30) 
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Proof. Using successively Proposition 13.31 (recall that 6 ~ h^), Lemma 15.21 and Theorem 13.81 we 
derive the first estimate: 

< Cth-'^\\I)N,(^i - c^iW-t + CoUi - ^\\o 

< CX(i-^)(||<Aillo + ll-^^ll*) < (5.31) 

To derive an estimate in the L°°(/2) norm, we proceed similarly: 

< (5.32) 
where we applied sequentially Proposition 13.31 Lemma 15.21 and Theorem 13.81 □ 

Proposition 5.9 For all t>2, there exists > such that 

l|Al,,,,.0i - ^L,,8<^h < Ct{\ loghf/' + h^-')\\cf>\U. (5.33) 

Further, for all e > and t >2, there exists Ce,t > such that 

l|Alg,,,,<Ai - Al,,,c/>^||oo,/. < Q,(/.*(i-«-(i+^)^ + h'-'mit. (5.34) 
Proof. We start by considering the decomposition 

In order to establish equations (I5.33P and (15.341) we first estimate the L°°{l2) norm (and, thus, 
the norm) for the first two terms of the right hand side of (j5.35p . For the first term, from 
Proposition 15.41 and (j5.ip we have 

< C\ \oghf/^{\\cpi - <p^\\o + \\Fn,<P^ - <P^\o) 

< a|log/i|3/2(/i*(i-/3) + /jt)||^||,. (5.36) 

For the second term, on the other hand. Proposition 15.51 and the fact that Gj(/>-' = show that for 
t > 2, 

IKAlg,,,, - Alg,,)Piv,0^||oo,/. < Cth<'-('^+^\\PN,<P'\\t + Ch^\\G,Pr,,<P'\Uj, 

< Cth'^'-''^+^\\<l>^\\t + Cih'-'\\cl^\\t, (5.37) 



32 



where in the last step we have applied the Sobolev embedding theorem H^^^ C L°°{l2) and (jS.ip . 
It thus remains to estimate the L°°{l2) and norms of the last term on the right hand side of 
equation (|5.35|) . 

We establish the bound first: in view of Proposition 13.31 and Theorem 13.81 we have 

W^L^M - n < C\W^ - < Cth'^'-^)U\\u (5.38) 

substituting this bound together with (|5.36p and (|5.37p in (|5.35p yields ()5.33p . To obtain the L°°{l2) 
bound we proceed similarly, but using the bound 

< Cr,th'^'-^^-^'+'^^\mt, (5.39) 

instead of (I5.38p . The first inequality in Equation (I5.39P follows from the Sobolev embedding 
theorem. In view of the relation 5 ~ , in turn, the second inequality results from (|3.24p with 
r := /3e/(l + /3) > 0. The third inequality, finally, follows from Theorem 13.81 and the fact that 
r + (1 + r)/3 = (1 + e)/3. The proof is now complete. □ 



6 A modified version of the Nystrom method 
6.1 Hermite interpolation 

As it has been presented and analyzed in the previous sections, the overall integral equation solver 
enjoys convergence of super-algebraic order for smooth data. The algorithm [5], however, incorpo- 
rates certain additional interpolation processes for added computational efficiency, as discussed in 
what follows. A convergence analysis for the modified algorithm is presented in Section [6.21 

One of the most computationally expensive segment of the Nystrom algorithm described in 
the previous sections concerns evaluation of the discrete operator Agj^^g ^ ^^^v^ : the high cost of 
this operation stems from the evaluation it requires of the trigonometric polynomial ^tVj at a large 
number of points 

{qhj,Z2 + {qhj - zi)tan0|) if |cos0^| < V2/2, (6.1) 

{zi + {q hj - Z2) cot 0^, q hj)^l~\ if | sin^^] < V2/2, (6.2) 

for all p = 0, . . . , 0j — 1 and all j = 1, . . . , J (see (|2.25p . ()2.27p and Figure [3|). The interpolation 
strategy mentioned at the beginning of this section accelerates this evaluation process. To describe 
the accelerated strategy in a manner that lends itself to analysis, we consider once again the operator 
L^-'^ in equation (|2.29p . we introduce the indicator function 

r 1, if Icos^l > V2/2, 
\ 0, otherwise, 

and we define L^^'j^ as the operator that is obtained as the angle-dependent weights c(0^) in (I2.29P 
are substituted by b{9i)c{6j,). Note that L^^/f contains only contributions from radial nodes lying 
on vertical grid lines (see Figure [3]); the remaining contributions to are collected in the operator 
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L'/h := ^Ih-'^lh- Accordingly, the operators A^^;^^^^ := {uj'or%l'fl and A^^;^^,^ := {uj'or%l'^ 
provide the decomposition 

^smg,S,h — ^smg,S,h ^ ^sing,<5,/i y^'"^) 

of the operator in equation (j3.13p . 

Each one of the operators on the right-hand side of ()6.3p lends itself to rapid (high-order) 
evaluation through application of univariate Hermite interpolation. To do this, given a positive 
integer d and a uniform mesh in the interval [0,1] with mesh-size k = 1/M, we consider the 
Hermite interpolation operator, that is, the operator : C'^[0, 1] — )• {(7 G C'^[0, 1] : g\[pk,{p+i)k] S 
IF'2d+i Vp} (where P2d+i is the space of polynomials of degree at most 2d + 1), which satisfies 

{Rk,dff\pk) = f^^\pk), p = 0,...,M, r = 0,...,d. 

The error estimate for this kind of piecewise polynomial interpolation is well known: 

11/ - Hfc,,/|Uoo(o,i) < C,fc2'^+2||^(2d+2) 11^^^^^^^ V/ e C'^'+'iO, 1]. (6.4) 

The Hermite interpolation method for approximation of the operator L^-'^^^ (resp. ^^^'j^O for & 
given j proceeds by interpolation of the trigonometric polynomial ^ via an application of the oper- 
ator Hfc^rf, in the variable ui (resp. U2), on a refined one-dimensional mesh of mesh-size k = hj/m 
{m > 1): the approximated operators thus result as the needed values of ^ at points of the form 
(j6.2p (resp. (j6.ip ) are replaced by the corresponding values produced by the Hermite interpolant 
^hj/m^dii'iQhj) (resp. ^hj/m4^{(l^ji'))- Note that this interpolation process, which requires eval- 
uation of ^ and some of its partial derivatives at points {phj/m,qhj) (resp. {qhj,phj/m)), can be 
implemented efficiently by means of the Fast Fourier Transform; see [5]. 

We express the Hermite-interpolation-based approximation of the operators L^-'^^ + L^-'^^ by 
means of the operator 

^iZ,m,d ■■= i^' ° r^) (l&'(H.,/™,. ^ I) + Lf,^(I ^ H,^/^,,)) . 

The corresponding matrix of operators for given values of h, m and d is denoted by ^^hej-m' 
that, for notational simplicity, the dependence on the parameters h, m and d is suppressed in the 
notation -4.^'^^^^,^^^. The new Hermite-interpolation-based discrete set of linear equations is given by 

S5Mrrn4>hcrm ■= [k^ + QhA'^'^VhVh + Qh^^herm^/^) </'herm = Qh^ ■ (6.5) 

In what follows we present a modified version of the analysis of Sections [3] through \5\ which 
takes into account the additional approximations associated with the Hermite interpolation process 
described above. This treatment thus accounts for all details of the basic partition-of- unity /polar 
integration Nystrom algorithm [5], and it establishes, with complete mathematical rigour, its con- 
vergence and high-order accuracy. 

6.2 Modified analysis, acounting for Hermite interpolation 

Analysis of the approximation (unique solvability, stability, convergence) for ()6.5p is carried out by 
comparison with the discrete solution cf)^. As we will show, m will have to grow as h decreases in 
order to retain convergence properties of the method. 
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Proposition 6.1 Assume that 6 ^ with 13 G (0,1) and that m = m{h) is chosen so that 
^-(2d+2)/^/3-i ^ 0. Then 

lim 11^5^ - Bs^herraWnO^nO = 0- 
/i->0 

Therefore, for h small enough the discrete equations ()6.5p are uniquely solvable and the operators 
hemi have uniformly hounded inverses. 

Proof. Note that the operators L^-'^^^ and L^^^^^ are weighted trapezoidal rules applied to 
V'm.isCP) + where V'u.iS is given by ()4.6p . Since iIju.s satisfies conditions (|4.4p . it 

foTiows that for all ^ e C°(72), 

Wi^t'Moo,!, + l|Lf;feiloo,/. < csuUj,. (6.6) 

By ()6.6p and the error bound for the Hermite interpolation ()6.4p we can estimate 

< ||Lf;^^(C - (H,^./„,, I)e)l|oo,7, + IlLff (e - (I » H,^./„,,)0||oo,7, 

< Q5/.f+^m-(2'^+2) (11^2^+2^11^^^^ + ^Clloo,/.). (6.7) 

Therefore, using the inverse inequalities ()A.23P and ()A.24p in the right-hand side of (j6.7p it follows 
that 

This inequality and the stability estimate for the interpolation operator Q/, (Lemma 15. ip imply 
that 

which proves the first assertion of the statement. Invertibility of Bs^herm as well as the uniform 
bound for their inverses follow from Theorem 13.71 □ 

In the last results of this section the constants will depend on /3 and d, without this dependence 
being made explicit. 

Proposition 6.2 Assume that 6 ^ h^ with (3 G (0, 1) and that m = m{h) is chosen so that 

^-(2a!+2)^/3-l _^ Q_ 2,et cj) and (/)herm ^6 the respective solutions of Bs4> = U and -B^^herm^herm = U. 

Then for every e > and t > s > with t > 1, there exists a constant Cs^t,e > such that 

U - 0hcrm||. < Cs,t,e{h'^^-^^-'\\ct>\\t + h"'^^^^-' m^^^''^^'' UU^^^,) . (6.8) 
Consequently if m{h) is chosen so that < Ch^"^^"^ for some r > 0, 

U - ^hermllo < Ch^''^'^+' \mi^2d+Z+r) / {l^p) ■ (6.9) 

Proof. Using (j6.7p . Lemma |5. II and the Sobolev embedding theorem (to bound (2d + 2)-th order 
partial derivatives by the f/'^d+s+c ^grm), we can easily prove that 

\\{BsMrm - Bs,hm\o < Ce5h'"'^^m-^^''+^)Uhd+^+,. (6.10) 
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Proceeding as in the proof of Theorem 13.81 (see (j5.25p and the bounds thereafter), we can introduce 
the discrete operator using Proposition 16.11 and bound 



ll*?!* — <?!'herm||o < C'||^(5,herm(<?!' " 0herm) ||o 

< C{\\V - QhUWo + \\(.Bs,h - ^5)0||o + ||(i35,herm " BsMo) ■ 

For the first two terms we use the same estimates as in the proof of Theorem 13.81 whereas the third 
term is bounded using (j6.10p . This gives (j6.8p with s = 0. For larger values of s we apply the same 
argument as in Theorem 13.81 (use inverse inequalities and the properties of Vh)- The estimate (|6.9p 
is a straightforward consequence of (16. 8p and the choices made for m. □ 



Theorem 6.3 Assume that 5 ^ with /3 G (0, 1) and that m = m{h) is chosen so that m (^'^+2) < 
(Jfi^-l^+r Jqj- some r > 0. Let tp be the solution of (|2.2p and V'herm be given by 

iGX(r) i=l j=l 

where 0herm = ('/'horrn' • • • ' '^herm)^'^ solution of Bs,herm4>hcTm = U. Then, denoting 

t := (2d + 3 + r)/(l-/3), 

we have the error estimate 

U - i^hermWmS) < Cr{\ log /l| § + /l'^^) || V'IIh^CS) ■ (6.11) 

Finally, for all e > 

H - V'hermllL-(S) < Cr,s{h'^'-^^-^'+'^^ + /l*"^) UWhHS)- (6-12) 

Proof. The proof of this result follows the same arguments as the one of Theorem 12.11 □ 

A Error analysis of the polar coordinate quadrature rules 

The discussion in this Appendix, which lies at the heart of the main convergence proof presented 
in this paper, provides estimates on the quadrature error 



oo f2n 



7p{p,6)dpd6 - Qh,k,-ti^ 



(A.l) 



-oo JO 

(which are used in Section 15. 2p for the family of polar-integration quadrature rules 

©-1 r oo -| „oo />27r 

Qm,7^:=A:^ c{ep)h ^{i{Gp) + q<Gp)h,ep) « / / i,{p,e)dpde (a.2 



p=0 



q= — oo 



with integrands arising from certain trigonometric polynomials. In detail, c{0) := min{l/| cos 0|, 
1/1 sin^l} (see ()2.28p ) and 7 is a certain piecewise smooth function which is typically discontinuous; 
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throughout this appendix we can and do assume that 7 G L°°(M) is a 27r-periodic function, without 
specific reference to its detailed smoothness properties. The integrand ip, in turn, is a function with 
compact support in the variable p which is given hy tp = xsS.N where, for a given trigonometric 
polynomial G Tjv, £,n{p,(^) '■= Civ(pcos 9, psin 0). 

The parameters in this quadrature formula are the positive integers and and the mesh- 
sizes h = 1/N and k = 27r/0; the angular and radial quadrature nodes, in turn, are given by 
Bp = pk ior p = 0, . . . , Q — 1 and 'y{Op) + qc{9p)h for q = —00, . . . , 00, respectively. Note that (|A.2p 
consists of trapezoidal rules in the variables and p, where the grid points used for integrating in 
p depend on 9. The error estimate (|A.1|1 is applied in Section [5] on the j-th coordinate patch for 
each j = 1, . . . , J — with N = Nj, Q = Qj and h = hj. 

Remark A.l In view of Proposition ^^ throughout this appendix the sequence {x5}s>o is assumed 
to satisfy the conditions listed in ()4.4p . Accordingly, the constants C and Cr in Theorem \A.l\ depend 
on the corresponding constants in the bounds (]4.4dp and on the constant cq in equation (|4.4cp . The 
dependence of the error bounds on the parameter 6, in turn, is made explicit in equations (IA.3P 
and (IXil) . 

The remainder of this section is devoted to the proof of the following theorem. 

Theorem A.l Let Q ~ A^i+" with a > and assume the sequence {x5}s>o satisfies the conditions 
listed in (14. 4p (see Remark \A.1\} . Then, there exists C > 0, independent of 5, h and 7 such that, 
for all 5 and h satisfying co6 < 1/2 and h < 6, we have 

£h,k,-y{xsCN)<c(h6-^\logh\ + 6^UN\\o VCiveTiv. (A.3) 
In addition, for each r > 2 there exists Cr > such that 

£h,kMsCN)<Crh'-6^-'-UN\\r V^^ G T^v- (A.4) 



A.l Error estimates for radial integration with trapezoidal rules 

This section provides bounds on the error 



£{^n;9) :-- 



ixsCN){p,9)dp-c{B)h iX5CN)i7{9) + qci9)h,i 



(A.5) 



that arises from the radial trapezoidal quadrature rule, where 7 is a given 27r-periodic function in 
L°°(M) and {xs} satisfies conditions (14. 4p . and where ^Ar(/5, 9) = iN{p§Jy9) ) for a given trigonometric 
polynomial G Ttv- Note that the mesh-size hc{9) and the location of the node corresponding 
to g = are allowed to depend on 0. As a first step we will limit our analysis to the case of 
trigonometric monomials, and we thus study the quantity 



Ern{9) := E(ern:^ 



{Xsern){p,9)dp-c{9)h J2 iX5enr){l{e) + qc{9)h, 



(A.6) 



where Smip, 9) = ^m{p§.{9)) = exp(27rijom • e{9)) 

For univariate functions defined on (—1/2, 1/2), we denote the Fourier coefficients by 



g{m) :-- 



1/2 



-1/2 



g{p) exp{—2iTimp) dp. 
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The Wiener algebra is the vector space of all functions whose Fourier series is absolutely (and 
therefore uniformly) convergent: 

oo 

n:={g:{-l,^)^C: ^ j^MI < oo}. 

m=—oo 

The Fourier series of g £ TZ gives a 1-periodic extension of the function to the entire real line. 
Lemma A. 2 Let f £ TZ be 1-periodically extended to the real line. Then the bound 



f{m) -h^f{-f + qh) exp(-27rim(7 + q h)) 

q=0 



<Y^\f{m + £N)\ 



holds for all m, 7 G M and N. 



Proof. Since the Fourier series of / converges absolute and uniformly, we can justify the following 
equality: 

N-i ^ N~l 

h ^ if exp{-2Trim-)){'y + qh) = ^ /(n) f /i ^ exp(27ri(n - m)(7 + g/i)) 

q=0 n£Z ^ cj=0 

= fin) exp(27ri(n — 771)7) ( ^ exp(27rig(n — m)h) 



However, 



and therefore 



N-1 00 

h (/exp(-27rim • ))(7 + qh) = ^ /(m + IN) exp(27ri£A^7), 

q=0 l=-oo 

which implies the result by separating the term corresponding to £ = in the right-hand side. □ 



Lemma A. 3 Let Em be given by (IA.6p . Then there exists a constant C such that 

Ern{0) < C{6 + h), V0 G [0, 27r] and Vm G Z^. (A.7) 
For all r >2 there exists Cr > such that 

:, if I tan^l > | cot 6\, 



, (N — \mi cot 0\ — \m2\Y ' 

Ernid)<CrS'-'- { 1 ^^'^ 

-— — j j j — — , otherwise, 

[JM — \mi\ — \m2tan0\y 

for all mi, m2 satisfying \mi\ + \m2\ < N — 1 and all 6 G [0, 2tt]. 
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Proof. Let Q denote the cardinality Q := ^{q : ^{6) + qc{9)h G {—cq5^cq5)}. Then, using ()4.4p 
and since, clearly, Q < f ^^jf ~'~ ^) ' have 



Ej,e) < \^^(p,e)\dp + cie)h\\xs\\Lo^iR2)Q<Co(2co5 + cie)h(^ + i 



-cqS 



c{e)h 



(A.9) 



and ()A.7p follows. 

To establish (|A.8p . in turn, we first consider the case 9 G [— 7r/4, 7r/4], so that c{9) = 1/ cos^ 
and I tan^l < 1 < | cot 9\. Clearly 



where 



/ 



c{9) 



exp(^27rimi-^^ dp-c{9)h ^ (/ exp(27rimi • )) + 9/^ 



q= — oo 



(A.IO) 



/('^) ■= Xsic{9)r,9) exp (27rim2c(6')r sin 0) = X5ic{9)r,9) exp (27rim2r tan 6*) 
Since co(5 < 1/2 and c{9) > 1, 

supvxs(.c{9)-,9) C (-co<5/c(e),co<5/cW) C (-1/2,1/2), 

it follows that the C°° function / is compactly supported in (—1/2,1/2), and therefore f £ TZ. 
Therefore, by Lemma lA.21 we have 



Ern{9) = C{9) 



q= — oo 



/(r) exp(27rimi r) dr — /i (/ exp(27rimi • )) + 



< ci9)Y,\K-mi+m)\. 

£^0 



(A.ll) 



Since for \mi\ + \m2\ < N — 1 we have | — mi + IN\ > \m2\ + 1 for all ^ ^ 0, it is possible to 
establish (jA.Sp from (jA.lip provided sufficiently restrictive bounds on |/(?7i)| for \m\ > \m2\ are 
obtained. 

To do this, letting b{r) := xs{i'-,d) and integrating by parts r times, we obtain 



/(m) 



b{c{9)p) exp(27ri(— m + m2 tan 9)p) dp 



-c{9) 



^ 2'K\{—m + 1712 tan 9) 
Using (14. 4p it therefore follows that 

l/(HI < ' ^ 



h^^\c{9)p) exp(27ri(— m + m2 tan 6)p) dp. 



\/27r 
1 



1 



m + m2 tan i 



\h^^\c{9)p)\dp 



1 



1 



y^vr y I — m + m2 tan 9^ c{9) J_ 



Co 5 



Co 5 



|6M(p)|dp 



< 



rl— r 



m + m2 tan 6*1 
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Inserting this inequality on the right-hand side of (jA.lip we obtain 

Em{0) < CrS^'^ y . ^— — -— . (A.12) 

— ^ ^ - ^\mi + m2tane - IN]"- ^ ' 

Since | tan0| < 1, the desired inequahty (lA.Sp follows from (1A.12P and Lemma lA. 8 1 This completes 
the proof for the case 6 G [— 7r/4, 7r/4]. For other values of 9 the proof follows similarly. □ 

Lemma A. 4 For all r > 2 there exists Cr > such that, for all mi, m2 satisfying \mi\ + \m2\ < 
N/2 we have 

Ern{0) < Cr5^-''h\ V0 G [0, 2^]. 
Proof. From Lemma lA.31 we obtain 

Em{e)<Cr6^- \^^ I I I ^ V0G[O,2^], |mi| + |m2| <iV-L 

— [N — \mi\ — \m2\Y 

The proof is completed by noting that for \mi\ + \m2\ < N/2 we have [N — \mi\ — \m2\)~^ < 
{N/2)-'- = Th''. □ 

In what follows we consider the 7r/2-periodic continuous even functions 

d{e) := min{|tan6l|,|cot6l|}, 

and 

^NiO):=[ iN-n-mdie))-^y^\ (A.13) 

(n,m)GH+ 

where 

H+ := {(n, m) : < n,m < N/2 <n + m < N}. (A. 14) 

Lemma A. 5 There exists C2 > such that the quantity E defined in (lA.SP satisfies 

£{Cn;0) < C2{s~'h + 6^'^N{e) + 6)UN\\o, V^^ G T^v, V0G[O,27r], 
Proof. Consider the mutually disjoint sets of indices 

Atv := |m = ?n'2) S : |mi| + |m2| < A^/2|, 

Hat := |?2], = (mi, 7712) G : |mi|, |m2| < A^/2 < |mi| + |m2| < 

= jin G : (|mi|, |m2|) G 
An := |m G : mi,m2 = ±iV/2|. 

We note that Ajv = for odd values of N and A^r contains four different elements for even N. 
Finally the set defined in equation (j2.32p satisfies Z|^ C A^v U Hjv U Ajv , with strict inclusion 
for even values of — since for even the set Z^ is not symmetric with respect to the coordinate 
axes in Z,"^ while all the index sets above are. 
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Given ^at G T^r we have 

£{im 0) < £i{^n; 0) + faC^Jv; 0) + £^{iN; 0), (A.15) 

where 

with the understanding that the last sum vanishes when A^r = 0. 

Using Lemma I A. 41 with r = 2 and the fact that the cardinahty ^A^v of A^v satisfies ^A^r < 
CiV^ = C/i"2 we see that that 

£i{^n;0)< meo^EmiO) V |6v(m)| < Cr^/iH^jvUo- (A.16) 

meA;v — 

From Lemma I A. 31 with r = 2, in turn, we have 

^ Ernief < C6-^ ^ (N - \mi\ - \m2di9)\r^ 

= Crt ^ (iV- |mi| - |m2(i(0)|)-^ (A.17) 

(To estabhsh the first inequality in ()A.17p use the transformation (mi, 7712) (m2,TO.i), which 
leaves Hjv invariant, in the cases for which d[6) = \ cot 9\.) Therefore 

£2{iN;0) < C6'^^N{0mN\\o. (A.18) 

Finally, for the remaining term we use Lemma [A . 3 1 and the fact that ^A^v < 4, and we thus obtain 

£3iCN;0)<C{5 + h)(^ ^ \CN{rn)\^y^^ <C{6 + h)UNh. (A.19) 
The lemma now follows from (|AT5il together with (jATTB]) . ijATTS]) and (|AT9]) . □ 



Lemma A. 6 For all r >2 there exists Cr > such that 

£{In]0) < Crh''6^-''\\^N\\r, S Tjy, V0 € [0,27r]. (A.20) 

Proof. The proof of this result proceeds through consideration of equation (lA.lSh . In view of 
Lemma I A. 41 we obtain 

£i{Cn;0) < Crd^~'-h'' \^N{m)\ < C',6^~'h''UN\\2, 

since for all ^ G we clearly have 

E ^11^112' C' = l+ n 12 M |2^2 • (A-21) 
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From (jA.lSP together with the fact that ^n{6) < 2 (as estabhshed in Lemma lA. 101 below) and the 
inequality 

1 /2 

I^MP) ^ /i'^(2\/2)nieilr y^GH"^ with r>0, (A.22) 

'mGZ2\A]v 

on the other hand, we obtain 

1 /2 

S2{^N;0)<Cr6~^(^ \^N{rn)\^) < C',6-^h'^UN\\r. 

Finally, using the first inequality in ()A.19p and (|A.22p we see that 

£3i^N;0)<C6(^ Y l6v(m)|')'^'<C'5/inie7v||r. 

This finishes the proof of the lemma. □ 



A. 2 Technical lemmas 

Lemma A. 7 In the space of bivariate trigonometric polynomials Tjv the following inverse inequal- 
ities hold: 

UN\\t<{V2hy-'UN\\s VCiveT^, s<t, (A.23) 

/ ^ \l/2 

lleivlloo,/^ < Yl \^N{ind\<N( Y l^^(l^)n =^"'ll^^llo VC^GT^v- (A.24) 
Proof. The proof follows directly from (j3.17p . □ 



Lemma A. 8 Let m,n G 7j be such that \m\ + \n\ < N — 1 and let r > 2. Then 

I < ^ o<e<i. 

\m + - iNl'^ ~ (N - \m\ - \n\CY ~ 



Proof. Noting that A^>A^— |m| — |n^| >A^— |n| — |m| >lwe obtain the chain of inequalities 
^ 1 ^ 2 2^ 1 

2^ \m nf: — fN\r — ( N — \m\ — \n\ ~^ 2^ 



\m + - tNl"- - {N - \m\ - \n\S,Y {IN - \m\ - \n\y 



CO 

< 



(N- 


\m\ — 






2 




{N- 


m — 






2 




{N- 


m — 





=2 

oo 

< 



1 + 



and the result follows. □ 
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Lemma A. 9 Let be given by ()A.14p . Then 

(n,m)e^+ 

Proof. Using the change of variables {j, m) = {n + m, m) the summation limits become 
(iV + l)/2 < j < iV - 1, j- N/2 <m< N/2. 

We thus have 



(N - j + mil - i)Y 



< 



E ( E 

-T^<j<N-l j-f<m<^ 



1 



< 



< 



V 1 N-j+l 

2 — ' — 

where the second inequality was obtained by using the fact that #{m : j — N/2 < m < N/2} < 
N — j + 1 and where the bound that introduces the integral relies on the fact that, for convex 
functions, the midpoint rule underestimates the value of the integral. 
To estimate Rn{€) we proceed as follows: 

dx 1 /l,^ ^ ^=1-2^ 



2 



< 



A2(l-e)2' 

This bound is vahd even when ^ = since Rn{'^) = 8{N - 1)/N^ < 8/N'^. 
Finally, noting that $Ar(i^) is an increasing function of ^ G [0, 1], we obtain 

r^-h dr r^-l dr 

and the result follows. □ 



Lemma A. 10 Let ^ n be given by ()A.13p . Then, there exists C > independent of N and @ such 
that 

e-i 

'^n{0)<2 and k'^'i'Nipk) < C{k + hlog{l/k)), 
where h = 1/N and k = 27r/0. 
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Proof. Note that = V^N ° d, which easily gives the first inequahty as a consequence of Lemma 
IA.9I A simple computation shows that for every integer i, 



I - d{e) > -\e - ce\ ee 

vr 



ce--,ce + j 



Ci 



(2^ + l)7r 



(A.25) 



Consider the sets 



IT vr 



\ [ci -k,Ci + k]. 



Then, letting 6p = pk, it follows that 



e-i 



p=0 



since the interval [q — fc, q + k] contains at most three of the grid points 9p and ^'jv < 2. Using 
Lemma lA.91 and ()A.25P we obtain 



V8 



By considering all possible arrangements of the uniform grid {Op : p = 0, . . . , — 1} with respect 
to the double interval Ii^k, it is easy to check that for all £ 



1 1 2 



E 



e 



< - ( 1 + lof 



® + <^(i + ioge). 



The proof now results by combining the three previous inequalities 

4^/27r 



e-i 

k^-^N{Op) <2Ak + 

p=0 



N 



:i + ioge) 



and collecting constants. 



□ 



A.3 Proof of Theorem [AJ] 



The first part of the proof does not rely upon the relationships between N and Q assumed in 
Theorem IA.lt up to a point in the proof that is clearly marked, N and are arbitrary positive 
integers, and, thus, k = 27r/0 and h = 2/N are unrelated real constants. The proof does assume 
throughout that h < 5, however. 

Let (/> G C~(M2) satisfy (9 + 27r) = (/>(p, 0) for ah p, 6* G R, as weh as supp (/>(-, 6*) C (-a, a) 
for all E M. A classical argument on the error of the trapezoidal rule for smooth periodic functions 
(see for instance [16^ Theorem 9.26]) shows that for all positive integers r we have 



< 



oo . fin 
-oo ^ Vo 



0-1 



0(p,^)de-A;^<A(p,ep))dp 

p=0 



dp < Ca^n 1 1 1 ioo(iR2). 



(A.26) 
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Noting that 

^KknixsiN) < Soixs^N) + kY^ £{^n; Op) (A.27) 

p=0 

and that (f) = X&£,N satisfies the hypotheses required to apply (|A.26P with a = cqS, the proof 
proceeds by producing bounds on the two terms on the right-hand side of ()A.26p . 

To obtain a bound for the second term on the right-hand side of (|A.27p we first invoke Lem- 
mas E5] and [XTTO] to establish that 

e-i _ e-i 
kY,£{iN]Op) < c{6-^h + 6-^k^^N{ep) + 6)UN\\o 

< C'(6-^h\og{l/k) + 6-^k + 6)UN\\o, Ve^GT^, (A.28) 



and we invoke Lemma IA.6I to show that for r > 2 we have 

e-i _ 

kY,£{^N;0p)<Ch'S^-'~UN\\r, y^N^TN. (A.29) 

p=0 

To bound the first term on the right-hand side of (|A.27p . on the other hand, we first note from 
the inverse inequalities ()A.23P and (|A.24p that, for all positive integers i, we have 

l|f^e^Af||L°°(R2) < ||(?a^iv||oo,/2 < Ceh^^ ^ ||9aCiv||o 

\a\<e \a\<e 

< c',h-^CN\\i<c';h-^-'UN\\o- 

(The first inequality follows from an application of the chain rule and the definition of S^n in terms 
of ^AT.) Therefore 



\\dg{X5^N)\\L'=°{R-^) < C^\\d^X5\\L°°{R^)\\d0 ^Cn\\l°°{R2) 

e=o 

r 

< C'UN\\oY.6-'h'-'-'<C"h-'~^UMh, (A.30) 

(the hypothesis 6^^h < 1 was used in the last inequality). Additionally, using ^AM\i and d'^h < 1 
we obtain 

H{XS^N)\\L^iR2) < Ch-'S-'^UNWr. (A.31) 

Inserting these inequalities into (jA.26p and since a = cq6 we obtain 

r 6h-^~-UM\\o, 

£oU5^N)<CkU (A.32) 

I 6'-^h''UN\\r. 

Using now the hypothesis G « N^~^°' with a > 0, it follows that k < Ch^^°' < Ch, and therefore, 
for r > 1/a, < Cah^~^^- Inserting the latter estimate in the inequalities (|A.32p we obtain 

SoiXs'^N) <CrV6^-''-UN\\r V^^ G T^y, (A.33) 
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for r = and for all r > 1/a, respectively. Therefore, equation ()A.33P holds for all r > by 
interpolation [19j. 

Using (jX27l), ()X28|) and the r = instance of (IXSgjl . and noting that \og{l/k) < Ci + 

C2log(l//i) equation (1X3]) follows. Equation ([Xil) . finally, follows from ([ATM]) and (Ia:291) (for 
r > 2). The proof is now complete. 
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